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ABSTRACT 



The Collapsar model provides a theoretical framework for the well known 
association between long gamma-ray bursts (GRBs) and collapsing massive stars. 
A bipolar relativistic jet, launched at the core of a collapsing star, drills its way 
through the stellar envelope and breaks out of the surface before producing the 
observed gamma-rays. While a wealth of observations associate GRBs with the 
death of massive stars, as yet there is no direct evidence for the Collapsar model 
itself. Here we show that a distinct signature of the Collapsar model is the 
appearance of a plateau in the duration distribution of the prompt GRB emission 
at times much shorter than the typical breakout time of the jet. This plateau 
is evident in the data of all three major satellites. These findings provide an 
evidence that directly supports the Collapsar model. Additionally, it suggests the 
existence of a large population of choked (failed) GRBs and that the 2 s duration 
commonly used to separate Collapsars and non-Collapasars is inconsistent with 
the duration distributions of Swift and Fermi GRBs and only holds for BATSE 
GRBs. 



Introduction 



There is a long line of e yidence connecting long GRBs (LGRBs) to col lapsing; massive 
stars (for recent reviewes see IWoosley fc Hegerll2006t iHjorth fc Bloomll201ll ). Among them 
are the association of half a dozens GRBs with spectroscopically confirmed broad-line Ic 
supernovae (SNe), and the identification of "red bumps" in the afterglows of about two 
dozens more, which shows a photometric evidence of underlying SNe. On top of that, the 
identification of LGRB host galaxies as highly star forming galaxies and th e localization 
of the LGRBs in the most active star forming regions within those gal axies ( iBloom et al. 
20021 iLe Floc'h etaDboO.j Ichristensen et all 12004 : iFruchter et aDl2006l ). provides an indi- 
rect evidence for the connection of LGRBs with massive stars. The model that provides 
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the theoretical framework of the LGRB-SN associati on is known as the Collapsar model 
(IMacFadyen fc Woosleylll999l ; iMacFadyen et al.ll200ll ). According to this model, following 
the core collapse of a massive star, a bipolar jet is launched at the center of the star. The 
jet drills through the stellar envelope and breaks out of the surface before producing the 
observed gamma-rays. However, although this model is supported indirectly by the LGRB- 
SN association, to this date we could not identify a clear direct observational imprint of the 
jet-envelope interaction, thus there is no direct confirmation yet of the Collapsar model. 

In this letter we analyze the expected duration distribution of the prompt GRB emission 
from Collapsars. We show that under very general conditions the time that the jet spends 
drilling through the star leads to a plateau in the duration distribution at times much shorter 
than the breakout time of the jet. We examine the duration distribution of all major GRB 
satellites and find an extended plateau in all of them over a duration range that is expected 
for reasonable jet-star parameters. We interpret these plateaus as the first identified imprint 
of the Collapsar model. Under this interpretation these findings (i) supports the hypothesis 
of compact stellar progenitors, (ii) imply the existence of a large population of chocked jets 
that fail to break out of the progenitor star and (iii) enable us to determine the transition 
duration that statistically separates Collapsars from non-Collapsars bursts, and to show 
that this time is individual for each satellite. The s eparation betwee n Collap sars from non- 
Collapsars is quantified and discussed in length in iBromberg et al.l (120 lid ). The letter is 
built accordingly: In section 2 we briefly review the propagation of a jet inside a star, and 
the duration of the prompt GRB emission. In section 3 we analyze the expected duration 
distribution of Collapsars. We compare the expected distribution with the observed one in 
section 4, and discuss our findings in section 5. 



2. The propagation of a jet in the stellar envelope 

As the jet propagates in the stellar envelope, it pushes the matter surrounding it, cre- 
ating a "head" of shocked matter at its front. For typical luminosities the head remains sub 
relativistic across most of the star (IMatznerl 120031 ; IZhang et al.l 120031 ; iMorsony et al.l 12007 ; 
Mizuta &: Aloyll2009l ; IBromberg et al.l 1201 lbl ). even thou gh the jet is ejected at relativistic 
velocities. The jet breaks out of the stellar surface after (IBromberg et al.l 1201 lal ): 

1/3 / Q \ 2/3 / n \ 2/3 / w \ 1/3 



tb — 15 sec • 



9 



R* 



15M m 



10 51 erg/sec y 

where Li SO is the isotropic equivalent jet luminosity, 9 is the jet half opening angle and we 
have used typical values for a long GRB. and M* are the radius and the mass of the 
progenitor star, where we normalize their value according to the typical radius and mass 
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inferred from observations of the few supernovae (SNe) associated with long GRBs. While 
the jet propagates inside the star the head dissipates most of its energy and the engine must 
continuously power the jet in order to support its propagation. For a successful jet breakout, 
the engine working time, t e , must be larger than the breakout time, t&. If t e < tb the jet fails 
to escape and a regular long GRB is not observed. 

Once the jet breaks out from the star it produces the o bserved emission at large distances 



from the stellar surface. Because of relativistic effects (ISari fc Piranl 119971 ) the observed 
duration of the prompt 7-rays emission, t y , reflects the time that the engine operates after 
the jet breaks out0: 

— t e — tb- (2) 

The distribution of i 7 is therefore a convolution of the distributions of the engine activity 
time and the breakout time combined with cosmological redshift effects. 



3. The duration distribution of Collapsars 

Under very general conditions, Eq. [2] results in a flat distribution of t 7 for durations 
significantly shorter than the typical breakout time. To see it, consider first a single value 
of tb and ignore, for simplicity, cosmological redshift and detector threshold effects. The 
probability that a GRB has a duration t 7 equals, in this case, to the probability that the 
engine work time is t 7 + tb- Namely, 

^(t-y)^ = P e (tb + tr^dtry, (3) 

where p 7 is the probability distribution of observed durations and p e is the probability 
distribution of engine working times. Expanding p e (t b + t 7 ) around t e = tb gives: 

Pe{tb + t 1 )=Pe(t b ) + 0(j£j (4) 

Thus p e (tb + ty) ~ p e {tb) = const for any t 7 that is sufficiently smaller than tb- Moreover, 
if p e (t e ) is a smooth function that does not vary rapidly in the vicinity of t e ~ tb over a 
duration of the scale of tb, then the constant distribution is extended up to times t 7 < tb- In 
the case of interest tb and t e are determined by different regions of the star: the breakout 
time is set by the density and radius of the stellar envelope at radii > 10 10 cm, while t e is 



1 Clearly the engine cannot be working in the same mode for a time that is longer than t 7 + % while if 
tb it is expected to l eave a clear signature on the temporal evolution of the prompt emission, which 



try t(i 

is not seen (Lazar et al 



200 



9) 
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determined by t he stellar core p roperties at radii < 10 8 cm. The core and the envelope are 



weakly coupled ( ICrowtherl 120071 ) and the engine is unaware whether the jet has broken out 
or not. Therefore, it is reasonable to expect p e (t e ) to be smooth in the vicinity of t e ~ £& 
and p y (tj) ~ const for t 7 < t b . In the opposite limit, where t b <C t 7 , then t 7 « i e , and eq. [3] 
reads: 

^ « I P f 6) , ^ , (5) 

In reality we expect t& to vary from one burst to another due to a variety of progenitor 
sizes and masses and a scatter in the jet properties. Moreover, effects such as the depen- 
dence of the break time on the jet luminosity may introduce correlations between t e and tb- 
Therefore we consider next a distribution of breakout times, Pb{tb), which can generally be 
correlated with p e {t e ) (we still ignore cosmological redshift and detector threshold effects). 
Since the stellar and jet properties are bounded (e.g., there is a minimal progenitor size and 
maximal jet luminosity), it implies that within the entire population of bursts there is a 
minimal breakout time, t^ m i n - Then eq. [3] is modified to: 

p 7 (t 7 )dt 7 = dt 7 / Pb{tb)Pe(tb + t^tbidtb , (6) 



where p e (t\tb) is the density distribution of t e given that the breakout time is tb- In this 
case all the arguments presented above for a single breakout time hold for each value of tb 
independently. Therefore, each distribution p e (t\tb) has a plateau at times t < t b with a 
normalization Pb(tb)p e (t b \t b )dt b . It can readily be seen that for t 7 < t b>min a plateau exists in 
all the distributions of engine work times, implying that the right hand side of Eq. [6] can be 
written as an integration over plateaus, with different normalization: 

Pl(tj ^ t b) min) ~ / Pb(tb)Pe(tb\tb)dt b = COnst. (7) 

Thus, p 7 has a plateau at times t 7 < tb,min- For any reasonable distribution of pb there is 
a typical value tf,, which dominates the contribution to the integral. Namely, bursts with 
t 7 < t b m i n are dominated by events with t b ~ t b . When p 7 is a monotonically decreasing 
function, as seen in the case of GRBs, then bursts with t b ~ tb dominate p 7 for any t 7 < tb 
and the plateau extends up to t b . 

Finally we take into considerations also the effects of cosmological redshift and detector 
threshold. All the quantities that we consider from here on are observed quantities. Namely, 
for bursts at a redshift z, all durations are stretched by a factor of (1+z), furthermore 
all distributions are the observed ones and are affected by the detector thresholds. For 
abbreviation we use the same notations as above. Let the observed redshift distribution be 
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p z (z)dz, which may be correlated with any other observed quantity. The observed burst 
duration distribution is then: 



p 7 (fc, 



p z (z)p b (t b \z)p e (t b + tJt b , z)dz dt b 



where p e (t\t b ,z) is the density distribution of t e given a breakout time t b and redshift z. 
Following the same considerations as above, if p e (t b \t b ,z) doesn't vary much in the vicinity 
of t b (for every t b and z), then p 7 has a plateau. Similar reasoning to the one that follows Eq. 
[7J implies that if p 7 is monotonically decreasing function, then a fiat distribution is observed 
up to t b , which is the breakout time whose contribution to the integral in [8] dominates its 
value at t 7 < t b . The integral is also expected to be dominated by bursts from a typical 
redshift z, implying that the typical intrinsic break time is t b /(l + z). 

At long observed durations, t 7 > t b , the distribution p 1 can in general depend on all 
three observed probability distributions, p b , p e and p z and their coupling. However, if bursts 
with the same breakout time, t b ~ t b , and at the same redshift, z ~ z, dominate the observed 
distribution both at t 7 < t b and at t 7 > t b , then p 7 (i 7 ^> t b ) oc p e {t^\t b ,z). Therefore, an 
extrapolation of p 7 (t 7 ^> t b ) to durations shorter than t b is similar to an extrapolation of 
Pe(te 3> tb\t b ,z) to a duration t e < t b . Namely, it is an estimate of the number of choked 
bursts. Now, if bursts with t b and z do not dominate the observed distribution at t 7 > t b 
then, for any reasonable conditions, the monotonically decreasing p 7 (t 7 > t b ) is decreasing 
less rapidly with t 7 than p e (t e 3> t b \t b , z). In that case extrapolation of p 7 to short durations 
(< t b ) again provides a reasonable estimate for the minimal number of choked GRBs. 

To conclude, under very general conditions the jet-envelope interaction in the Collapsar 
model is predicted to produce a plateau in the duration distribution of GRBs at short 
observed durations. This is true for any breakout time, engine working time and redshift 
distributions (including cases where the various distributions are correlated) as long as the 
engine working time distribution is smooth enough. The bursts in the constant section of 
p 7 are dominated by a population with an observed breakout time t b , and the plateau is 
extended up to i 7 ~ t b . Based on the typical observed GRB parameters (see eq. [T]) and 
given that a typical GRB is observed at redshift k 2 we expect t b m 50 s. Finally, it is 
well established that at much shorter durations (< 1 s) the duration distribution contains 
a significant population of bursts that are not ass ociated with the death o f massive stars, 



which are known as short hard GRBs (SGRBs) ( Kouveliotou et al.l Il993l ; iNarayan et al. 



200 ll ; lMatznerll2003l : iFox et al.ll2005l : iBerger et al 



2005 



Nakarl 120071 ). These bursts are non- 



Collapsars, and the above arguments don't apply to them. Therefore, when considering the 
overall burst duration distribution we expect a flat section for durations significantly lower 
than 50 sec down to the duration where these non-Collapsars dominate. 
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Fig. 1.— The T 90 distribution, dN/dT 90 , of BATSE (red), Swift (blue) and Fermi GBM 
(green) GRBs. Also plotted is the distribution of the soft (hardness ratio < 2.6) BATSE 
bursts (magenta). For clarity the Swift values are divided by a factor of 5 and the Fermi 
GBM by 15. The dotted line that ranges down to m 2 sec mark the duration range where 
Collapsars constitute more than 50% of the total number of BATSE GRBs. At shorter times 
the sample is dominated by non-Collapsars. Note that the quantity dN/dT is d epicted and 
not dN/d\ogT as traditionally shown in such plots (e.g.. lKouveliotou et al.lll993l ). The black 
lines show the best fitted flat interval in each data set: 5 — 25 sec (BATSE), 0.7 — 21 sec 
(Swift), and 2.5 — 31 sec (Fermi). The upper limits of this range indicate a typical breakout 
time of a few dozens seconds, in agreement with the prediction of the Collapsar model. The 
distribution at times > 100 sec can be fitted as a power law with an index —4 < a < —3. 
Soft BATSE bursts show a considerably longer plateau (0.4 — 25 sec), indicating that most 
of the soft short bursts are in fact Collapsars. 



-7- 



4. The observed distribution the prompt GRB durations 

The observed duration of a GRB is characterized using T90 ~ t 7 , during which 90% of 
the fluence is accumulated. We use the data from the three major GRB detectors: BATSE, 
Swift and Fermi GBM. For BATSE we use the current catalog (04/21/91-05/26/00; contain- 
ing 2041 bursts). The data of Swift is taken from its online archive] (12/17/04-08/27/11: 
containing 582 bursts). Fermi data is extracted from GCNs using the GRBox websitqj 
(08/12/08-07/21/11; containing 194 bursts). Each data set is binned into eq ually spaced 



logar ithmic bins, where the minimal number of events per bin is limited to five (IPress et al. 



19891 ) . A bin with less than five events is merged with its neighbor. We use a x minimiza- 
tion to look for the longest logarithmic time interval that is consistent with a flat line within 
1 a, where the only free parameter is the normalization. We verify that varying the bin size 
doesn't change the length of the plateau by much. 

Fig. 1 depicts the observed distribution of T 90 , p 7 (T 90 ), for the three major GRB 
satellites. Note that we show here t he quantity jP 7 (Tgn) = dN/dT and not dN/dlogT 



traditionally shown in such plots (e.g.. iKouveliotou et al.lll993l ). The best fitted flat regions 



are highlighted in a solid bold line on top of each distribution. In all satellites these plateaus 
range about an order of magnitude in durations (BATSE 5-25 sec, 3.6/6 x 2 /dof; Swift 0.7- 
21 sec, 8.9/7 x 2 /dof; Fermi 1.2-31 sec, 4.1/6 x 2 /dof). The extent of the plateau varies 
slightly from one detector to another. This is expected given the different detection threshold 
sensitivities in different energy windows (see below). At the high end of the plateau the T90 
distribution decreases rapidly and can be fitted at long durations (>100 s) by a power law 
with an index, a, in the range —4 < a < —3. 

The existence of the plateaus and their duration range (~ 2 — 25 s) agrees well with 
the expectation of the Collapsar model. However, one cannot exclude the possibility that 
the origin of the observed flat sections is unrelated to the effect of the jet breakout time 
that we discuss above. For example, these plateaus may somehow arise coincidentally from 
the combination of two distributions: one increasing (LGRBs) and one decreasing (SGRBs). 
In order to test the hypothesis that the observed plateaus indeed reflect the distribution of 



LGR Bs, which are Collapsars, we use the fact that SGRBs are harder (IKouveliotou et al. 



19931 ). Restricting the analysis only to soft bursts should preferentially remove SGRBs from 
the sample. Thereby, LGRBs should dominate the duration distribution of a sample of soft 
bursts down to durations that are shorter than in the case of the entire sample. Now, if 



2 http://swift. gsfc.nasa.gov/docs/swift/archive/grb_table 
3 http://lyra. berkeley.edu/grbox/grbox.php 
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LGRBs are collapsars and there duration distribution is flat at short times, then T90 distri- 
bution of a sample of soft bursts should exhibit a plateau that extends to shorter durations 
than the T go distribution of the whole sample. We present in fig. 1 also a distribution of 
BATSE soft bursts (magenta), which are defined as bursts with hardness raticQ, HR < 2.6, 
the median value of bursts with T90 > 5 sec. Remarkably, the plateau in this sample ex- 
tends from 25 sec down to 0.4 sec (15.2/12 % 2 /dof), over almost two orders of magnitude 
in duration, compared to the original 5 sec in the complete BATSE sample. This lands a 
strong support to the conclusion that the observed flat distribution is indeed indicating on 
the Collapsar origin of the population. It also implies that HR is a good indicator that 
effectively filters out a large number of non-Collapsars from the GRB sample. 



5. Discussion 

The observed plateaus in all three durtion distributions, and most notably in the distri- 
bution of the soft Batse bursts, provide a direct support for the Collapsars model for LGRBs. 
An inspection of different regions of the observed temporal distribution (Fig. 1), under the 
interpretation of the plateau as an imprint of the time it takes the jet to break out of the 
envelope, provides further important information. 

1. The end of the plateau and the decrease in the number of GRBs at long durations, 
allows us to estimate the typical time it takes a jet to breakout of the progenitor's en- 
velope. All three distributions are flat below ~ 10 sec in the GRB frame, implying that 
tb ~ a few dozen seconds. This value fits nicely with the canonical GRB parameters 
taken in Eq. [lj an d provides anot her support that the stellar progenitors of Collapsars 



must be compact (jMatznerl 120031 ) . 



2. The end of the plateau at long durations is dominated by the distribution of jet break- 
out times rather than by the engine working times distribution. Thus, as we discuss in 
section [31 an extrapolation of the distribution at long durations, i 7 3> tb, to durations 
shorter than tb provides a rough estimate of the number of Collapsars with engines 
that don't work long enough for their jets to break out. If the jet fails to break out 
then it dissipates its energy into the stellar envelope producing a so called a "choked" 
or "failed" GRB. There have been numerous suggestions that such a hypothetical pop- 
ulation of hidden "choked GRBs" exists and that these are strong sources of high 



4 HR is the fluence ratio between BATSE channels 3 (100-300 keV) and channel 2 (50-100 keV). 
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gravitational waves (INorris 



2003 



energy neutrinos (lEichler fc Levinsonlll999l; iMeszaros fe Waxmanll200ll ) and possibly 



Daigne fc Mochkovitchl 120071 ) 



At long durations p 7 (T 90 > 100 sec) can be fitted well with a power law, p 7 oc Tg with 
—4 < a < —3. Extrapolating to t e < t^ we find that if p e continues with this power 
law to t e < tb there are many more chocked GRBs than long ones. For example, even 
if we extrapolate this distribution only down to t e = t&/2 there are still ~ 10 times 
chocked GRBs than long GRBs. This prediction is consistent with the suggestion that 
shock breakout from these choked GRB s produces a low lum inosity smooth and soft 
GRB jTan et al.lboOll ; IWang et~aTll2007t Eakar fc Sarilboilh . Indeed the rate of such 
low l uminosity GRBs is far larger than that of regular long GRBs ( jSoderberg et al 



20061 ). 



Eichler et al. 


1989; 


Fox et al. 


2005; 


Berger et al. 


2005; 


Nakar 


2007) 



classify since all their hard energy properties largely overlap with those of the Collapsars 



(jNakarll20071 ). The least overlap is in the duration distributions and hence, traditionally 
a burst is classified as a non-Collapsar if T 90 < 2 sec. Even though this criterion is 
based on the duration distribution of BATSE it is widely used for bursts detected by 
all satellites. 

The fraction of Collapsars at short durations can be estimated by extrapolating the 
plateau in the duration distribution. Since there is also an overlap with the SGRBs at 
long durations, the real hight of the plateaus are somewhat l ower than what i s show n 
in fig. 1. A detailed treatment of this issue is presented in iBromberg et al.l ( 1201 lcT ) . 
Nevertheless we can use the observed plateaus to obtain a crude estimation to the 
duration below which the majority of GRBs are non-collapsars, by extrapolating the 
plateaus to the bins in which Collapsars constitutes 50% of the bursts. This shows that 
for BATSE the transition occurs at ~ 3 s, while for Swift it occurs at ~ 0.7 s and at 
~ 1.2 s for Fermi, although the statistics in the Fermi data is quite poor. The shift in 
the transition time is expected since non-Collapsar bursts are also harder and different 
detectors have different energy detection windows. BATSE has the hardest detection 
window, making it relatively more sensitive to non-Collapsar GRBs. Swift has the 
softest detection window making it relatively more sensitive to Collapsar GRBs. 

Thus, putting the dividing line between Collapsars and non-Collapsars at 2 sec is 
statistically reasonable for BATSE bursts. However, our analysis shows that it is 
clearly wrong to do so for Swift, and likely so to Fermi bursts. 



4. While the difference in the lower limit of the flat ranges is understood qualitatively 
in view of the detection windows of the different detectors, the variance in the upper 
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limit is less obvious. It may reflect various selection effects in triggering algorithms. 
Unfortunately Swift's complicated triggering algorithm makes is difficult to explore 
this point quantitatively. A more interesting possibility is that it reflects a physical 
origin, e.g. that different satellites probing populations with different This could 
be explored when the statistical sample of Fermi GBM becomes sufficiently large. 

To conclude we remark that it is intriguing that the unique feature discovered here for 
the GRB's duration distribution depends just on the simple relation (Eq. [2]). This implies 
that a similar plateau is expected in any transient source whose duration is determined by 
two unrelated processes in a similar manner. In particular, in any case in which one process 
emits radiation and another independent one blocks it for a while. For example, this would 
take place in a source engulfed by dust which will be obscured until the dust is destroyed 
by the radiation wave. This could be of importance in interpreting numerous transient 
observations, particularly now at the dawn of astronomy in the time domain. 

This research was supported by an ERC advanced research grant, by the Israeli center 
for Excellent for High Energy AstroPhysics (TP), by an ERC and IRG grants, and a Packard, 
Guggenheim and Radcliffe fellowships (RS), by the ISF grant No. 174/08 (EN) 
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